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We introduce a finite volume renormalization scheme for the ,/V-Majorana-component O(A0 in- 
variant Gross-Neveu model. Universal observables are defined that are accessible to precise nu- 
merical simulation in various discretizations and allow for an extrapolation to the continuum limit. 
Here first numerical results with Wilson fermions are reported. For N — 2 they reproduce exact 
finite volume continuum results in the massless Thirring model. Our ,/V = 8 data are ready for 
comparison for instance with staggered results in the future. 
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1. Introduction 

In lattice QCD several fermion discretizations are in use in current dynamical simulations. 
Beside the very costly recent variants with chiral symmetry, Wilson and staggered fermions are 
the standard choices with their well-known relative merits and weaknesses. For the latter choice, 
unphysical multiples of 4 degenerate flavors (in 4 dimensions) can only be avoided by the 'rooting' 
procedure which is under much debate [R1], as universality of the continuum limit is not guaranteed 
any more by the locality of the action. 

We thought that in this situation a two-dimensional study, where the continuum extrapolation 
can be controlled better than in QCD, would be a valuable check. Gross-Neveu models (GN) [Q] 
and the Schwinger-model come to mind. The latter is closer to QCD in being a gauge theory, while 
the renormalization structure of GN is more realistic as its coupling is dimensionless, even asymp- 
totically free for N >3, instead of being superrenormalizable. For us the latter aspect prevails. 

The work of the ALPHA collaboration has demonstrated, that a precise continuum extrapo- 
lation becomes feasible for quantities like the Schrodinger functional of QCD, where the system 
size is used as a physical scale to probe the field theory. We hence construct a similar finite vol- 
ume renormalization scheme for GN. Also here the finite size supplies an infrared scale allowing 
the mass to be tuned to a critical value 1 corresponding to the (here discrete) chirally symmetric 
continuum limit. With the coupling as the only remaining free parameter, the situation becomes 
similar to the massless QCD Schrodinger functional. While the aim clearly is to use this scheme 
also for staggered fermions, possibly with 'rooting', in the future, at present we only have Wilson 
simulations to report on. 



2. Model and renormalization scheme 

We consider the action density of the Euclidean theory 

Here the Grassmann spinor £ caries a two-valued spin and an Af-valued flavor index such that 
an internal O(A0 symmetry arises and the antisymmetric matrix ^ obeys ^y^^ 7-1 = — /X. With 
this symmetry no other 4-fermion interaction is possible and the model is renormalizable as it 
stands. This is in contrast to the chiral GN model with N Dirac fields, which shares continuous 
chiral symmetry with QCD, but allows for two independent couplings and a third one with Wilson 
fermions [g]. This makes it much more difficult to approach a definite continuum theory and we 
hence restrict ourselves for our fermion-testbed to the 0(/V) invariant class. For m = the action 



( |2.1| ) has the additional discrete invariance 

^75^ (2-2) 

that breaks spontaneously for N = °° [§]. 



It vanishes if the chiral symmetry is preserved by the regularization. 
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Our action (2.1) has so far referred to the continuum. In the Wilson r = 1 lattice regularizaton 
we have to replace 

dfi-^dn, m^m--d ll d* (2.3) 



with the forward (d), backward (d*) and symmetric (d) lattice derivative. Now ( PL2| ) is not a 
symmetry at finite cutoff. It emerges however in the continuum limit for a suitable tuning m = 

m c (g 2 ). 

For our finite volume scheme we consider the field £ on a T x L torus with (anti)periodic 
boundary conditions in space (time) 

£(jc + r6) = -§(*), £(*+Ll) = +£(*), (2.4) 

i. e. a spatial ring at finite temperature. In the following we take T = L, aspect ratio one, and the 
smallest momentum then is p* = {%/T,Q). In a momentum version of our finite size scheme we 
formulate a complete set of normalization conditions on 2- and 4-point functions using external 
momenta ±/>*. We however here prefer to use correlations at physical separations in space-time 
for numerical reasons and also in the hope - supported by perturbation theory - to minimize cutoff 
effects in this way. 

We Fourier-transform in space only 

|(*0,/>)=a£e-^(*), !(*,) = Ifo.O) (2.5) 

and impose normalizations on the 2-point function 

0= (| T (r/4)^|(0)V (2.6) 



is a renormalized field, and the first condition is required by ( |2.2D and determines 
the critical mass m c . Note that at separation T /2 the analogous equation holds as a consequence 
of antiperiodicity and time-reflection invariance for all m and could not serve to define m c , while 
with our choice we found good sensitivity to do so. On the lattice only T/a that are multiples of 
4 must be simulated to obtain a scaling situation, an acceptable restriction. Finally a renormalized 
coupling is obtained from 

£r = ^1 (Ui(T/2)n R A0) g 2 (r/2)^| S)2 (0)) , (2.8) 

where the subscripts on E,r refer to two specific flavor values. For Wilson fermions we have per- 
formed a 1-loop calculation and obtain 

am c = -(N - l)Kg 2 + 0(a 4 g 2 ,g 4 ), ^ = 0.384900179460 (2.9) 
Z^ = l + 0(/) (2.10) 

g 4 + 0(g 6 )) (2.11) 



^ = ^(| T (r/2)^ 7 ol(0)). (2.7) 



8l = ¥ ^rj8 2 + 



(N _ 2){ M^ +co]+Ci+OUl] 
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with 

c = -0.483524477, c y = 0.30965176. (2.12) 

The 1-loop value for m c agrees with the literature [|], ^j. The coefficient of the logarithm in 
C p- 11 ) has the well known value and vanishes for the Thirring case N = 2. Note that a finite 
coupling renormalization as well as the linearly divergent mass renormalization are still there. 



3. Simulation 

The standard approach for GN is to factorize the 4-fermion term with a scalar auxiliary field 
e £«*E,«W = f Yld^G)^ 2 ^^. (3.1) 

J x 

Usually a Gaussian field is employed giving 2 

Z= fl\dae-^ al [PfA(a)] N , (3.2) 

J X 

where the Pfaffian results from integrating out £ . The operator 

A(o)=^{y ll d tl +m + gO -ad*d) (3.3) 

can be taken real antisymmetric in the Majorana representation of y„. Hence for even N we have 
the non-negative weight < [PfA] w = [detA T A] A '/ 4 which may be represented by N/2 real pseud- 
ofermions for an HMC approach. 



4. Massless Thirring model (N = 2) 

We have extended the well known exact continuum solution of the Thirring model to our 
finite geometry and plan to report on this elsewhere This allows to predict many correlation 
functions. As usual, correlations of the U(l) current 

Jn{xo,Pi)=J dx^-^yy^x) (4.1) 

are particularly easy to obtain, where we use one Dirac field here at A 7 = 2. In this case the sym- 
metry (^2|) gets automatically promoted to an axial U(l). Correlations depend on the continuum 
coupling g cont . By varying it we produce the curve in Fig. [jj Note that the quantity on the ;t-axis is 
the thermal expectation value of the squared total charge on our ring. 

In the lattice transcription we employ the exact Noether U(l) current that does not renormalize. 
We take the continuum limit by extrapolating from L/a = 16 . . . 48 for values g = 0.4 and 0.7 of the 



bare lattice coupling. The mass is tuned to m c defined by (2.6) on each lattice. The extrapolations 
for g = 0.4, leading to the lower point in Fig. [j] are shown in Fig. ^[ The other point is similar but 
with larger errors. 
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Figure 1: Exact universal relation between current correlations in the massless Thirring model (N = 2). 
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Figure 2: Continuum extrapolations producing the lower point in Fig. [I]. 

5. Gross Neveu model (N = 8) 



Here the coupling renormalizes similarly to QCD. We now adjust the coupling ( J2.8| ) to the 
value gR = 0.38 and take the continuum limit as discussed before. Our universal observable in this 
case is 



(M«( x o,Pi 



1 



NZ, 



a^- ipiXl (e(0WU(x) 



(5.1) 



for scalable xq values and admissible p\. Two examples are shown Fig. [3] and Fig. |j. 
6. Remarks on staggered fermions 



For naive Majorana lattice fermions, i. e. leaving out the (oc a) part in ( |2.3| ), we find the usual 
taste multiplicity 2 D in D dimensions (each momentum component around or Tt/d). If one now 
tries to 'spin-diagonalize' by transforming £ one can only achieve a reduction factor 2 D / 2 ~ ! in 



"At finite N only a finite number of moments in rj are relevant allowing in principle also other distributions 
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Figure 3: Continuum limit of a Gross-Neveu correlation at N = ! 
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Figure 4: As in Fig. |J but for a different time separation. 

contrast to iPl 1 for Dirac fermions when yf, \[r are changed independently. The reason is that V 
can only be reduced to 2 x 2 blocks in the Majorana case but not diagonalized. Since the Majorana 
form is natural for GN we simply deal here with naive fermions and 4 tastes (on top of N flavors). 
The 4-fermion interaction term naively would read in 2-momentum space 

7^4 I d^^ Pi )i i T (p l )n i (P2)lj(P3)nj(P4)- 

PU ...,P4 " * ' 

taste— mixing 

As indicated, there would be contributions with for instance p\ and p% in different corners of 
the Brillouin zone which mix tastes in the bilinears that are flavor-scalar. This problem was al- 
ready solved in an early effort to simulate the model [Q]. An additional factor Yl tl cos(a(p\ + 
P2)u,/2) cos(a(p3 + p4)n/2) under the sum is expected to enforce taste symmetry in the continuum 
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limit. In position space this corresponds to distributing the interaction term over a plaquette. In 
this form one then naively expects an additional taste symmetry in the continuum limit which may 
reproduce results with exact flavor symmetry in the Wilson formulation. 

Acknowledgement: We would like to thank Bjorn Leder and Peter Weisz for critically reading this 
manuscript. 
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